A Gl-like globular cluster in NGC 1023 ^ 

S0ren S. Larsen 
UC Observatories / Lick Observatory, University of California, Santa Cruz, CA 95064, USA 

soerenSucolick . org 

ABSTRACT 

^^ ' The structure of a very bright (My = —10.9) globular cluster in NGC 1023 is 

CN ' analyzed on two sets of images taken with the Hubble Space Telescope. Prom careful 



^ 



X 



1. Introduction 



modeling of King profile fits to the cluster image, a core radius of re = 0.55 ±0.1 pc, 
effective radius R^. = 3.7 ± 0.3 pc and a central surface brightness of fJ^oiV) = 12.9 ± 0.5 

rn , mag arcsec"^ are derived. This makes the cluster much more compact than lu Cen, 

but very similar to the brightest globular cluster in M31, Gl = Mayall II. The cluster 
^ I in NGC 1023 appears to be very highly flattened with an ellipticity of e « 0.37, even 

\^ I higher than for co Cen and Gl, and similar to the most flattened clusters in the Large 

O ' Magellanic Cloud. 

I> 

Subject headings: galaxies: star clusters — galaxies: individual (NGC 1023) 
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^ ' A very bright globular cluster in the nearby lenticular galaxy NGC 1023 was recently identified 

by Larsen, Brodie and Kissler-Patig (2001). At a distance of 9.9 ib 0.6 Mpc (Ciardullo, Jacoby and 
Harris 1991), NGC 1023 is close enough that images taken with the Hubble Space Telescope (HST) 
will show globular clusters (GCs) as spatially extended objects, and the structure of individual 
clusters with sufficient S/N can be studied in considerable detail. It is of interest to compare 
the structural parameters (size, ellipticity etc.) of globular clusters in different galaxies because 
significant differences are known to exist even among GCs in the Local Group. For example, clusters 
in the Large Magellanic Cloud (LMC) are more flattened than those in the Milky Way (Geisler & 
Hodge 1980). It has also been noted that the brightest GCs in the Milky Way and M31 (a; Cen 
and Gl = Mayall II, respectively), as well as in the Magellanic Clouds, are the most flattened in 
their respective host galaxies (van den Bergh 1984), and van den Bergh (1996) has suggested using 
the HST to test if the brightest globular clusters around other galaxies are also flattened. In the 
present paper the structure of the bright globular cluster in NGC 1023 is analyzed using archive 
HST images, and it is then compared with to Cen and Gl. 



^Based on observations obtained with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope 
Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc. under NASA 
contract No. NAS5-26555. 



2. The bright cluster nl023-13 

Larsen, Brodie and Kissler-Patig (2001) presented spectroscopy of 11 old globular clusters in 
NGC 1023, most of which were selected from WFPC2 exposures obtained as part of Program 6554 
(PI: Brodie) and discussed in Larsen and Brodie (2000). However, a few additional objects were 
selected from ground-based images, one of which (nl023-13) turned out to be a very bright globular 
cluster. The 2000.0 coordinates are a = 2:40:27.84 and 5 = 39:04:40.2. Ahhough this object is 
close to the saturation limit on the Brodie dataset and located near the edge of one of the WF 
chips, it is also included on a shorter WFPC2 exposure of NGC 1023, obtained for D. Richstone 
(HST Program ID 6587). This dataset consists of 5 integrations of 260 sec each in the F555W band 
and 2 integrations of 900 sec in the F814W band, short enough that the cluster is well below the 
saturation limit. 

Fig. 1 shows a close-up of the Richstone WFPC2 image of the cluster, compared to a star of 
similar brightness. The cluster is clearly resolved and appears to be significantly elongated. At 
a distance of 9.9 Mpc, one WFPC2 pixels spans 4.8 pc, comparable to or somewhat larger than 
the typical half-light radius of a globular cluster (e.g. Kundu and Whitmore 2001), so special care 
must be taken when analyzing the light profile of clusters at this distance. Nevertheless, the S/N 
of nl023-13 is sufficient to allow a relatively detailed analysis of its structure and compare with 
other bright globular clusters like lo Centauri in the Milky Way and Gl = Mayall II in M31. 



2.1. Modeling of the light profile 

Of the 5 available F555W exposures in the Richstone dataset, only two were used for the 
analysis presented here. The two selected exposures required no shifts in the x and y directions, 
and consequently no resampling of the Point Spread Function (PSF), before combination. Addition 
of the individual images, including elimination of cosmic ray hits, was done with the IMCOMBINE 
task in IRAF^. Inspection of the data quality files for the Brodie WFPC2 image showed that the 
image of nl023-13 was, in fact, not saturated, so we are in the fortunate situation of having two 
WFPC2 images of the same cluster, allowing for a useful consistency check of the results. The 
cluster is located on two different WF chips (WF2 on the Richstone and WF4 on the Brodie 
dataset), and at very different positions within the chips {{x,y) = 210,560 versus {x,y) = 740,79) 
in the two pointings. 

The F555W images of the cluster were modeled using the i shape algorithm (Larsen 1999), 
assuming various analytical models for the intrinsic luminosity profile of the cluster and then 
convolving with the HST PSF, generated by the TinyTim software (Krist and Hook 1997). A 
convolution with the WFPC2 "diffusion kernel" was also included. The ellipticity, Full Width At 



^IRAF is distributed by the National Optical Astronomical Observatories, which are operated by the Association 
of Universities for Research in Astronomy, Inc. under contract with the National Science Foundation 



Half Maximum (FWHM) and orientation of the model profiles were iteratively adjusted until the 
best match with the observed image was obtained. As the TinyTim manual and WFPC2 handbook 
give two slightly different diffusion kernels, fits were done for both of them. Because the diffusion 
kernel is wavelength dependent but only given for the F555W filter in the manuals, no fits to the 
F814W images were attempted. 

Two different types of analytical profiles were assumed: King (1962) profiles of the form 

2 

, r <rt (1) 



/x(r) = k 
and "Moffat" profiles given as 
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The King profiles are characterized by a core radius Tc and a tidal radius rt and are known to 
provide excellent fits to globular cluster luminosity profiles. The Moffat profiles are similar to the 
profiles used by Elson, Fall, & Freeman (1987) to fit young LMC clusters, but note that the a in 
Eq. (2) is equivalent to their 7/2. In the limiting case of a = 1, Eq. (2) is identical to a King profile 
with infinite tidal radius. 

The version of ishape used here allows the concentration parameter c = rt/rc of the King 
profiles and the exponent a of the Moffat profiles to be kept at a fixed value, or to vary these pa- 
rameters during the fitting process. The c parameter is the most uncertain of the fitted parameters, 
and the exact value returned by ishape was quite dependent on the initial guesses. Generally, the 
fitted c values were between 100 and 300. In order to ensure better stability in the fits and assess 
the effect of variations in c, the remaining shape parameters (ellipticity, FWHM, orientation) were 
thus determined for fixed c values of 100, 200 and 300. The Moffat fits, on the other hand, were 
quite stable and consistently returned a values around 1.20 it 0.05. 

The fitting radius was set to 15 pixels or l'.'5, corresponding to a radius of 72 pc for a distance 
of 9.9 Mpc. A number of fits were also done for a fitting radius of 25 pixels, but for this larger radius 
the background gradient becomes noticeable, although the results of the fits remained essentially 
identical to those done with the smaller radius. 

Figure 2 shows the results of various fits to the Richstone image. From the left: King profiles 
with c = 10,30,100,300 and cx) and the best-fitting Moffat profile with a = 1.20. The top row 
shows the best-fitting model images generated by convolution of the King/Moffat profiles with 
the TinyTim PSF and diffusion kernel, and the bottom row shows the residuals when the model 
images are subtracted from the observed image. The King profiles with c < 100 or c = 00 are clearly 
unable to reproduce the observed cluster image, while the c = 100 and c = 300 profiles show little 
structure in the residuals beyond the photon noise. Formally, the best fit (lowest x^) '^ss obtained 
for c ~ 200. Increased S/N could potentially help in constraining the King profile c parameter and 
for this purpose a number of fits were carried out on all 5 Richstone images combined. These tended 
to converge towards c ~ 150. It thus seems reasonable to adopt c = 200 it 100 for the concentration 



parameter. The a = 1.2 Moffat profile also provides a very good fit, but the King profile fits which 
are physically better motivated will be used for most of the discussion in this paper. 

The various King and Moffat profiles are plotted in Fig. 3, scaled to the same luminosity 
within I'.'S and computed for the actual fitted core radii. Note that Fig. 3 shows the intrinsic 
cluster profiles, not the observed ones (which are obtained by convolving the profiles in Fig. 3 with 
the HST PSF and diffusion kernel). Since 1 WF pixel spans (/.'I, it is clear that the central parts of 
the profiles are unresolved and the differences in the fits rely mainly on the behaviour in the wings. 



2.2. Shape parameters for best fitting King profiles 

The best fit parameters are given in Table 1 for each of the two WFPC2 datasets. Each value 
in the Table is an average of 6 individual fits, for King profiles with c = 100, 200 and 300, and for 
the two different diffusion kernels. The quoted errors are simply the standard deviation of these 
individual measurements, which is probably a better estimate of the true uncertainty than the 
standard error on the mean. 

The first rows in the Table give the photometry in a 0^5 (5 pixels) aperture, corrected for 
reddening towards NGC 1023. The photometry was done with the PHOT task in DAOPHOT 
directly on the images, measuring the background as the mode of all pixels in an annulus starting 
at 20 pixels and 10 pixels wide. Calibration to standard V and / magnitudes was done following 
Holtzman et al. (1995) and using the Schlegel et al. (1998) extinction value of Ab = 0.262 mag. 
The distance modulus of 29.97 it 0.14 derived by Ciardullo, Jacoby and Harris (1991) was based on 
Burstein and Heiles (1984) extinctions, which often tend to be somewhat lower than the Schlegel 
et al. (1998) values. However, for NGC 1023 the two happen to agree to within 0.01 mag in Ab- 

For an extended object like nl023-13, the O'.'S aperture slightly underestimates the total bright- 
ness. According to Holtzman et al. (1995), the WFPC2 PSF itself scatters about 10% of the hght 
to radii beyond their 0'.'5 reference aperture, while measurements of nl023-13 through an r = 20 
pixels aperture gave total V magnitudes brighter by about 0.2 mag compared to the r = 5 pixels 
measurements. Thus, the V magnitudes in Table 1 are probably too faint by ~ 0.1 mag, and after 
applying this correction the absolute magnitude becomes My = —10.9. Likewise, direct integration 
of the King profiles also yields a correction of about 0.1 mag from r = 0'.'5 to r = oo. 

There is excellent agreement between the photometry from the two pointings, as well as between 
the ishape profile fits. The two sets of fits give FWHM values of 0'.'027 ± 0'.'006 and 0'.'029 ± O'.'OOT, 
minor/major axis ratios of 0.62 zt 0.02 and 0.64 it 0.01, and the position angles on the sky agree 
to within 3 degrees. Note that the relatively large range in the c parameter does not lead to 
correspondingly large uncertainties on the fitted FWHM. 

The FWHM and concentration parameter returned by ishape can be converted to more famil- 
iar quantities such as the core radius {vc) and half-light (effective) radius {Re)- For a King profile 



the FWHM and r^ are related as 



FWHM = 2 




(3) 



i.e. FWHM ~ 2rc for c ^ 1. A similar simple analytical relation between the FWHM or re and 
Re does not exist, but can be approximated by a power-law of the form 

i?e/rc~ 0.547 cO-4S6 (4) 

For c > 4 this approximation is good to ±2%. 

These relations assume that the profiles are circularly symmetric, while the FWHM values 
returned by i shape are measured along the major axis of the fitted profiles. In order to compute 
the core and effective radii in Table 1, an average of the major and minor axis FWHM has been 
used. This results in core radii of Vc = 0.52 it 0.12 pc and re = 0.57 ± 0.13 pc for the two datasets, 
and effective radii of Re = 3.47 it 0.33 pc and Re = 3.85 it 0.34 pc, respectively, at the assumed 
distance of NGC 1023. 



2.3. Ellipticity 

The minor/major axis ratio returned by the fits is quite stable at around 0.63 it 0.01, corre- 
sponding to a very large ellipticity of e = 0.37 ± 0.01. An independent estimate can be obtained by 
using the ELLIPSE task in the STSDAS package to fit elliptical isophotes directly to the cluster 
image, although these will tend to underestimate the ellipticity because the true shape is blurred 
by the HST PSF. This effect will be worse closer to the center, while ellipse fits at large radii are 
uncertain because of low S/N. At r = 5 pixels, the ELLIPSE task returns ellipticities of e = 0.24 
on both the Richstone F555W and F814W exposures, and 0.14 and 0.20 on the Brodie F555W 
and F814W exposures. Irregularities in the background due to the proximity to the edge of the 
image are clearly visible in the Brodie exposures and the ELLIPSE fits are probably less accurate 
for these data. 

ELLIPSE fits were also done on the model images generated by ishape (Figure 2). These fits 
returned an ellipticity of e = 0.25 at r = 5 pixels, in nearly perfect agreement with the fits to the 
actual science images. Thus, it is clear that the blurring of the cluster image by the HST PSF leads 
to an underestimation of e when measured directly on the images, and the true answer probably 
lies closer to the e ~ 0.37 value from the King profile fits. 



2.4. Central surface brightness 

Although the central parts of the King profiles are unresolved on the WFPC2 images, the 
central surface brightness for the cluster can still be computed, assuming that the King profile 



continues all the way to the center. The central surface brightness hq of a King (1962) profile and 
the total luminosity L{R) within a radius R < rt are related as 



and 
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For reddening-corrected V = 19.19 and V = 19.17 within a radius of i? = Cf.'b, this results in 
/io = 12.8 lb 0.4 mag arcsec"^ and fiQ = 13.0 it 0.4 mag arcsec"^ for the two fits. Again, the errors 
are computed as the standard deviation of the 6 values based on the individual fits. The computed 
/io values are not very sensitive to the exact concentration parameter and most of the uncertainty 
comes from the core radius. Alternatively, one can also compute the central surface brightness for 
the Moffat profile fit. This yields fiQ = 13.6 mag arcsec^^, somewhat fainter than for the King 
profile. However, a King profile is probably a better approximation to the true light profile of the 
cluster and a value around /xq = 12.9 it 0.5 mag arcsec"^ may be a reasonable estimate of the central 
surface brightness and its associated error. 



3. Comparison with Gl and uj Cen 

Table 1 also lists a number of properties for two other bright globular clusters: u Cen, the 
brightest GC in the Milky Way, and Gl = Mayall II, the brightest known GC in M31 (Mayah 
and Eggen 1953). Pritchet & van den Bergh (1984) first noted a significant flattening of Gl and 
estimated an ellipticity of e = 0.22. Using HST images of Gl, Rich et al. (1996) measured a core 
radius of re = O'.'lTOitO'.'Oll or 0.63 ±0.04 pc for a distance modulus of 24.42 (Preedman & Madore 
1990) and found a tidal radius of 28'.'2 = 105 pc, corresponding to a concentration parameter c = 166 
and an effective radius of 4.1 pc. They obtained a central surface brightness of 13.5 mag arcsec^^ 
in V and a slightly higher ellipticity (e = 0.25 it 0.02) than the one reported by Pritchet & van den 
Bergh (1984). Also using HST images, Meylan et al. (2001) measured a somewhat smaller core 
radius of re = 0.52 pc, but a larger tidal radius of 200 pc for Gl. They quote a half-mass radius 
of rfi = 14 pc, corresponding to an effective radius of 10 pc (Spitzer 1987, p. 16) if light traces 
mass, but their Table 3 indicates an effective radius closer to 1'.'18 or 4.4 pc, in better agreement 
with Equation (4) which gives Re = 0.52 pc x 0.547 (200/0.52)°-'^^'5 = 5.1 pc. Contrary to Rich et 
al. (1996), Meylan et al. (2001) found a variation in ellipticity with radius, ranging from e = 0.1 in 
the innermost and outermost parts, to e = 0.3 at a radius of 2^ or 8 pc, with a mean of e ~ 0.2. 
All of the above studies agree on an integrated V magnitude of My = —10.9. 

Most of the relevant parameters for to Cen are listed in the McMaster catalogue (Harris 1996). 
This cluster has My = -10.24 and is also quite elongated with e = 0.19 (Frenk k Fall 1982). It 



is, however, more extended than Gl and has a core radius of 3.8 pc. The effective radius is 6.2 pc, 
somewhat larger than the typical 3 — 4 pc for globular clusters. 

Compared to these two bright globular clusters in the Local Group, nl023-13 appears to be 
nearly an identical twin of Gl. Both clusters are of comparable luminosity and their effective radii 
are quite similar to those for globular clusters of much lower luminosities. This is consistent with 
the lack of correlation between effective radius and cluster luminosity observed in the Milky Way 
(Djorgovski & Meylan 1994). Both Gl and nl023-13 have very compact cores and, as a result, very 
high central surface brightnesses. In fact, the highest central surface brightness for any Galactic 
globular cluster is /Uq = 14.15 (NGC 1851, from the McMaster catalog) and uj Cen itself has a much 
lower ^0 = 16.8 mag arcsec"^. 

From the core radius and central surface brightness, the central density po can be estimated 
from the relation 



3.44 X 10^° 



--0-4W){V')/'A/f/r^M^T,n-3 



Po = —B lO-^-^''"^^nWi)M0pc-3 (7) 

Pre 

with P « 2 and re in pc (Peterson & King 1975; Wilhams & Bahcah 1979). For Gl and nl023-13 
this leads to po = 2.0 x 10^ Mq pc^^ and po = (4.1 ± 2.4) x 10^ Mq pc~^, respectively, assuming a 
M/L ratio of 1.6 (Illingworth 1976). Here we have used the average of the two estimates of po for 
nl023-13 listed in Table 1. From King-Michie models, Meylan et al. (2001) obtain a central density 
of Po = 4.7 X 10^ M0 pc~^ for Gl, in fair agreement with the estimate obtained from Eq. (7). 
Again, nl023-13 has a similar high central density to Gl, while lv Cen has a central density of only 
~ 1400 Mq pc~3^ 

The metallicity of uj Cen is given as [Fe/H] = —1.62 in the McMaster catalogue. Metallicity 
estimates for Gl range between —1.19 it 0.25 (Brodie and Huchra 1990) and —0.7 (Rich et al. 
1996), while Meylan et al. (2001) find an intermediate value of [Fe/H] = -0.95 ± 0.09. For nl023- 
13, Larsen, Brodie and Kissler-Patig (2001) found [Fe/H] = —1.15 it 0.22 from Keck spectroscopy, 
and a very similar value of [Fe/H] = —1.11 from the V — I color. nl023-13 thus seems to be 
somewhat more metal-rich than uj Cen, but of comparable metallicity to (or maybe slightly lower 
than) Gl. It is also worth noting that - unlike most globular clusters - both ui Cen and Gl have 
quite a large internal scatter in their metallicity distributions. Whether or not this is the case for 
nl023-13 is, however, impossible to tell from the current data. 

One common property of all these bright globular clusters is their high ellipticities. The ishape 
fits indicate an ellipticity as high as e ~ 0.37 for nl023-13. uj Cen has an ellipticity of e = 0.19, 
while Gl appears to display a range of ellipticities as a function of radius, reaching a maximum of 
e = 0.30 at intermediate radii. The present data do not permit an analysis of internal variations in 
the shape parameters (ellipticity, orientation) of nl023-13, but from the mean values it seems that 
the cluster may be even more fiattened than uj Cen and Gl. Is an ellipticity e = 0.37 unrealistically 
large? Although most globular clusters in the Milky Way have low ellipticities, clusters in the 
Large Magellanic Cloud tend to be much more elongated. One of the more striking examples is 
the intermediate-age (~ 2 Gyr) cluster NGC 1978 with an ellipticity of e = 0.30 (Fischer, Welch, 



& Mateo 1992). After correction for projection effects, Geisler &: Hodge (1980) find that very few 
LMC clusters are likely to be spherical, most having true ellipticities in the range 0.2 - 0.4. Thus, 
the high ellipticity for nl023-13 might not be unreasonable, especially if the cluster is seen nearly 
edge-on, but it would definitely rank it as one of the most elliptical clusters known. 

Which mechanism might be responsible for the high ellipticities of massive globular clusters? In 
a recent paper, Cen (2001) has argued that many globular clusters may initially possess significant 
amounts of angular momentum, which should cause flattening of their initial shapes (King 1961). 
Some evidence has been presented that rotation might indeed be a dominant mechanism in causing 
flattening of globular clusters (Davoust 1986; White & Shawl 1987) and of uo Cen in particular 
(Meylan & Mayor 1986; Merritt, Meylan, & Mayor 1997). On the other hand, Fischer, Welch, k 
Mateo (1992) found no evidence for rotation in the highly flattened cluster NGC 1978, although 
they could not rule it out. However, it is not quite clear why this effect should affect the most 
massive clusters more strongly. Perhaps the tendency for more massive clusters to have longer 
relaxation times (Spitzer 1987, p. 40) can cause an initial flattening to persist for a longer period 
of time in these systems, or maybe the more massive clusters tend to be born with larger angular 
momenta. Detailed dynamical calculations and a better understanding of globular cluster formation 
will be required before these questions can be definitively answered. 

Another mechanism which could produce elongated clusters with high angular momenta is 
merging of binary clusters (Sugimoto & Makino 1989). If the two clusters have different metallicities 
then the result might be an elongated cluster with a metallicity spread, as observed in ui Cen and 
Gl, although it seems unlikely that two clusters born near each other (in time and space) would 
exhibit such large metallicity differences. Alternatively, the merger might result from a later chance 
encounter between two clusters (Icke & Alcaino 1988), but this clearly makes the mechanism less 
attractive as a general explanation. 

Finally, Meylan et al. (2001) have discussed the possibility that Gl and uo Cen are the remnants 
of nucleated dwarf galaxies. This might make the metallicity spread and elongated morphology of 
these two clusters more understandable, but does not appear to be a likely explanation for the high 
ellipticities of most LMC clusters. 



3.1. Sources of errors 

Comparison of the two independent fits gives a good estimate of the random errors involved 
in the fitting process. However, because the characteristic sizes involved are on the same order of 
magnitude as the resolution of the images, the fitted parameters are strongly dependent on correct 
modeling of the WFPC2 PSF. One critical component is the "diffusion kernel", which takes into 
account charge diffusion between neighboring pixels in the CCDs. The diffusion kernel given in the 
TinyTim manual scatters less light to the surrounding pixels than the one in the WFPC2 handbook, 
and this difference also influences the fits. By comparison of the fits with the two different diffusion 



kernels, it was found that the kernel given in the WFPC2 handbook generally gave smaller FWHM 
values by about 0.02 - 0.04 pixels or 0'.'002 — 0'.'004. This is consistent with this kernel scattering 
more light into surrounding pixels, thus causing i shape to fit a slightly narrower intrinsic profile. 

We also fitted the star shown in Fig. 1 with and without the diffusion kernel. When using 
the diffusion kernel, ishape was unable to improve the King model fits relative to comparison fits 
obtained by using a delta function, as expected for an unresolved stellar image. In contrast, if the 
diffusion kernel was omitted then ishape clearly recognized the star as an "extended" object with 
an effective radius of about 0.2 pixels, confirming that the "raw" TinyTim PSF underrepresents 
the size of a point source and that the diffusion kernel needs to be taken into account. 

Errors in the distance to NGC 1023 are not included in the uncertainty estimates in Table 1. 
The 9.9 lb 0.6 Mpc value from Ciardullo, Jacoby and Harris (1991) is based on the planetary nebula 
luminosity function, and agrees well with the 10.2 Mpc distance estimate given by Faber et al. 
(1997). Size parameters (FWHM, Vc, Re and rj) scale linearly with the distance, while the central 
surface brightness /xq is distance-independent. Thus, from Eq. 7, the central density is inversely 
proportional to the distance. Compared to the measurement errors, the uncertainty on the distance 
to NGC 1023 has only a minor effect on the results, even if the ±6% error estimate by Ciardullo, 
Jacoby and Harris (1991) might be somewhat optimistic. 



4. Summary 

Based on a detailed analysis of its spatial structure, the bright globular cluster nl023-13 is 
found to be very similar to the Gl cluster in M31. Both have an absolute V band magnitude of 
My = —10.9. Using the average values from Table 1, the spatial profile of nl023-13 is well fitted 
by a King model with core radius 0.55 it 0.1 pc, effective radius 3.7 it 0.3 pc and concentration 
parameter c ~ 200. The cluster appears to be very flattened with an ellipticity around e ~ 0.37, 
making it one of the flattest star clusters known, and providing another example of the van den 
Bergh (1984) observation that the brightest GCs in galaxies also tend to be the most flattened. 
Assuming that the King profile continues all the way to the cluster center, a very high central 
surface brightness of lJ-o{V) = 12.9 it 0.5 mag arcsec^^ is derived, implying a central density of 
(4.1 lb 2.4) X 10^ M0 pc~^. This is comparable to Gl, but two orders of magnitude higher than for 
ijj Centauri. 

This work was supported by National Science Foundation grant number AST9900732. I thank 
Jean Brodie and Sidney van den Bergh for comments on the paper and enlightening discussions, 
and the referee for useful suggestions. 
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Fig. 1. — Images of the globular cluster nl023-13 (left) and a star (right). The two objects are of 
about the same integrated magnitude and are shown with the same contrast settings. The cluster 
image is clearly more extended and appears to be elongated. Each image spans 40 x 40 pixels or 
4" X 4". 




Fig. 2. — Synthetic cluster images and residuals for various analytical cluster models, fitted by 
ishape. From left to right: King (c = 10, 30, 100, 300, cxd) and Moffat {a = 1.20). 
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Fig. 3. — Illustration of the various luminosity profiles fitted by ishape. King models with con- 
centration parameters c = 10, 30, 100, 300 and oo (dotted/dashed lines) and a MOFFAT profile 
with index a = 1.20 (solid line). All profiles have been normalized to the same luminosity within 
r = I'.'S. Note that these are the intrinsic profiles, not taking into account the fact that the actual 
observed profiles are blurred by the WFPC2 PSF. 
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Table 1. Comparison of nl023-13, Gl and uj Cen. 







Richstone 


Brodie 


Gl" 


UJ Cen^ 


Vo 




19.190 ±0.005 


19.173 ±0.005 


_ 


_ 


{V - /)o 




1.038 ± 0.006 


1.050 ±0.006 


- 


- 


Mv 




-10.9 


-10.9 


-10.9 


-10.24 


[Fe/H] 




-1.15 ±0.2 


-1.2. ..-0.7 


-1.62 


c 




200 ± 100 


166. ..380 


17.4 


minor/major 


0.62 ±0.02 


0.64 ±0.01 


~0.75 


0.81 


P.A. 




-40.9 deg ±1.1 


-38.3 deg ±1.0 


- 


- 


FWHM 




0'.'027 ± 0'.'006 


0'.'029 ± 0'.'007 


- 


- 






1.27 ±0.29 pc 


1.39 ±0.32 pc 


- 


- 


Tc (pc) 




0.52 ±0.12 


0.57 ±0.13 


0.52. ..0.63 


3.8 


Re (PC)^ 




3.47 ±0.33 


3.85 ± 0.34 


4.5 ±0.5 


6.2 


n (pc) 




95 ±27 


105 ± 30 


90. . . 200 


67 


^0 {V mag 


arcsec ^) 


12.8 ±0.4 


13.0 ±0.4 


13.5 


16.8 


po (M© pc- 


-^) 


(4.7 ±2.8) X 10^ 


(3.5 ±2.0) X 10^ 


4.7 X 10^ 


1.4 X 10^ 



Note. — For r^ and Re the given values are averages of minor and major axis values. Vq and 
iy — I)q are given for a 0'.'5 (5 pixels) aperture. 

c = King profile concentration parameter = rt/vc. P.A. = Position angle on sky, measured N 
through E. yUo = central surface brightness, pQ = central density. 

"": Tabulated values for Gl represent the range quoted in the literature. See discussion in main 
text for details. 

: data from the McMaster catalog (Harris 1996) 
'^: Computed from Vc and rt- 



